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Abstract - The report deals with the possible 
implications of pressure diffusion for shocks in one 
dimensional traveling waves in an ideal gas. From 
this new hypothesis all aspects of such shocks can be 
calculated except shock thickness. Unlike conven- 
tional shock theory, the concept of entropy is not 
needed or used. Our analysis shows that temperature 
rises near a shock, which is of course an 

experimental fact; however, it also predicts that very 
close to a shock, density increases faster than 
pressure. In other words, a shock itself is cold. 


Introduction 

In inviscid flow theory, fluid acceleration arises from 

dv 

a gradient of pressure by p(~-g)--Vp ; here we reverse 

dv 

this idea and regard -p(~- g) instead, as a pseudo -pressure 

gradient which arises from acceleration. Thus in this 

unconventional approach, pressure is regarded as a special 

type of energy density, the gradients of which contribute to 

fluid accelerations. As a type of energy density, pressure 

might diffuse (not a violation of energy conservation) 

and we postulate that it does so at a rate determined by the 

dv 

difference between Vp and -p( ~-g) , that is, at a rate 

dv 

V* [ k { Vp+p (— - g ) } ] , where k is a coefficient of diffusion. 

For viscous fluids this diffusion term is generally non- 
zero. In the absence of accelerations, pressure gradients 


2 


would undergo Laplace diffusion. 

The implications of the concept of pressure diffusion 
for fluid dynamics seem quite general. Pressure diffusion 
might be a prominent factor in the irreversibility of 
compression or expansion of an ideal gas, a fundamental 
phenomenon which conventional fluid dynamics is unable to 
explain or predict from first principles without recourse to 
ad hoc extreme increases in a poorly understood thermodynamic 
term, the second coefficient of viscosity [1, pp . 304-308]. 

Furthermore, pressure diffusion is consistent with the 
experimental observation that the speed of sound in CO 2 , N 2 O 
and SO 2 at 1E5 Pa and 313° K increases with frequency [2], 

In fact the coefficient of pressure diffusion k (taken, like 
all other thermodynamic terms in this report, as a function 
of density p and pressure p) for air near standard conditions 
can be estimated from such experimental data [3] to be 
300 m^/sec . 

In this new approach to fluid flow modeling, the state 
variables of the model are density p, the three velocity 


components in the vector v, and pressure p. To the usual 


Navier- Stokes mass conservation equation for tt* and 

o t 

d v 

momentum conservation equation for T~ is added to new 

o t 

equation for . The pressure equation can be derived from 


an energy conservation equation in which pressure diffusion 
is present as a term and from specification of an internal 
energy function for the fluid. 
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Fundamental Equations 

The fundamental equations of Newtonian fluid dynamics 


with pressure diffusion are [1, pp . 2, 47-49] 


(i) 


and 


(v* V) v+p' 1 [ -Vp+ M V 2 v+(f+^)VD 
+(M+M*) Vp-jDVp+DVf ]+g 


( 2 ) 


Here p is dynamic viscosity; f is the second coefficient of 
viscosity; M is the matrix with {dv^SXj} in row i, column 
j; M denotes the transpose of M, the matrix with {av J /3x i } 
in row i, column j; and g is a constant gravitational 
acceleration vector. Obviously we require that p and f are 
differentiable functions of p and p. Dilation D is the trace 
of the matrix M, that is D-tr (M ) L . Implicit in (2) 
is the assumption that a special property of the energy 
density p is that the gradient Vp contributes to 

o t 

In this report a fluid is regarded as a continuum of 
matter which can be completely described by its density p 
velocity v and pressure p. A fluid is presumed to change 
with time in such a manner that matter and momentum are 
conserved; thus use will be made of two of the familiar 
Navier - S tokes equations (1) and (2). Density, velocity, and 
pressure dynamics can then be predicted provided an equation 
for ^7 in terms of density, velocity and pressure 
(including the spatial gradients thereof) can be obtained. 
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This is done in the next section. Indeed it is shown that 

[ [/xtr{M(M+M*) } + (r-f/OD 2 
+V« ( t VT ) +V • [ « { Vp+p (^r) ■ g ) ) ] ] (3) 

where T is temperature (a function of p and p), r is a heat 
diffusion coefficient, k is a pressure diffusion coefficient, 
and B is the bulk modulus of elasticity (all functions of p 
and p). Throughout we assume that p and p always have 
positive values and that the parameters u, p , f, r and k are 
analytic functions of p and p; these parameters are all non- 
negative and is assumed positive. If /c~0 , then (3) 
reduces to a well known but seldom used energy conservation 
equation [4, p. 33]. 

3. Derivation of the Pressure Equation 

The region R of flow will be assumed to be finite 
although possibly changing with time and for which at each 
time a general version of Gauss' Divergence theorem holds. 

The boundary dR of R must change smoothly with time and at 
each time 3R must be a closed oriented surface with an 
outward pointing normal vector field n defined uniquely 
almost everywhere. It follows that for any smooth vector 
field F over R, 

j V • F dV - I F • n dA . 

Jr idR 

Let us regard R as partitioned into a large number of small 
cubical cells, each with edge length i. 
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A standard expression of conservation of energy in 

Newtonian fluids [4, p. 33] is 

^* t (pu+^pv*v)--V* [ ( pu+p+hpv* v) v] + V« [ (T -j^ODv 

4*/i (M+M* ) v ] +V • (rVT)+pg*v (4) 

The term - V* [ ( pu+p+Hpv* v) v] arises from the rate of energy 

carried by fluid flow through the walls of a fixed cell, 

namely (pu+p+^pv»v)vn, where n is the outward normal 

vector of the cell. The fluid carries energy as internal 

energy per unit volume pu and kinetic energy density ^pv*v. 

The rate of work through cell walls by viscous forces results 

in the next term. The term V* (rV T) arises from the flow 

of energy through the cell walls by the conduction of heat, 

and p g*v is, of course, the contribution of gravitational 

force to the rate of change of energy in the cell. 

Pressure is the result of a force per unit area. In 

the presence of acceleration, a momentum change can be 

equated with a force. Thus if a fluid accelerates through 

the wall of a cubic cell of edge length i, the pseudo pressure 

generated by the momentum change (force) is pi 3 (^-g)/i 2 , 

including the effect of gravity. Discounting gravity, the 

pseudo pressure gradient (between two such cubes which share 

dv 

a face) generated by the momentum change is p>(^T-g). Thus 

the net rate of pressure diffusion in general form is taken 

dv 

to be V* [ k { Vp+p g) ) ] > where k is a diffusion parameter. 
This term is thus added to the list of power flux summands on 
the right-side of (4) to yield 
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d 2 * 

~( pu+hpvv)--V» [ ( pu+p+hpv . v) v] +V • [ (f-y/O Dv+p (M+M ) v ] 

+V* ( t VT ) + pg • v+V • [/c{Vp + p(yy-g) } ] (4a) 

combining (2) and (4a) leads to 

P^--pD+Mtr[M(M+M*) ] + (r-p)D 2 

+V. (rVT)+V. [/c{Vp + p(^--g) ) ] (5) 


s ince 


and 


du_iu djg. £u 
dt 3p dt * 3 p p 


^-^ +v .Vp 

dt at p 


Ae. 


we can use (5) to solve for 77 , that is, 

o t 

^■--V Vp+(|y) ' 1 [ pf^ - ’^] D+ (pf^) ' 1 [/*tr{M(M+M*) ) 
+ (f -|p)D z +V. (rVT)+V. [«{Vp+p(yy-g) } ] ] 


( 6 ) 


If we can show that the bulk modulus of elasticity B-p(~^) s , 
where s is entropy [4, pp . 165] can be written in terms of 

the internal energy u of the fluid by the relation 

then making use of this relation into equation ( 6 ), one 
obtains equation (3). This is shown as follows. 

Generally u-u(p,p). If specific entropy is held 
constant, then we may write p~f(p) and 

Recall that P“P 2( |^)s a consequence of the Gibbs equation 


du 


Tds+p~ 2 pdp. It follows that 


E. _ flu „ ,&E.x , 

p dp p ( ap )s p d P 


and so 




- O’ 1 


[-p 


3u 

dp 


(7) 
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as claimed. In fact (7) can be used to define B without 
reference to entropy. Thus the concept of entropy is not 
essential in describing the dynamics of Newtonian fluids. 


3. One-Dimensional Flow 

For one - dimens ional fluid flow, the equations (1), (2) 

and (3) reduce to 


$JL 

at 

£v 

at 




= - V 


ax 
3v 
ax 


(pv) 


+ p~ M 


flE. , 5 / 


ax 


(«)] + g 


ax'''ax 


All 


. b £v iu r 1 f 2 

at ax B ax (p ap ) [ ^ ( ax ) 

. 2_. IT ^_r ,dv _ . ni 

dx^ 3 x' 3 x^ dx ®) ] ] 


.Ae. 


»£v 


- v*— - B y r + ( a 1 f n ( ® V ) ^ 

ax ax ^ap ; l ^ax' 

, fl_, aXs ^ fl_ r flUs 1 1 

+ ax (, dx> + ax lK ax ( *ax )11 


-StLs 


( 8 ) 

(9) 

( 10 a) 

( 10 b) 


where rj - pu + f . 

For an ideal gas 


u 


ITY- so B=yp and p (|^) * 1 -y-l , 


P ( 7 - 1 ) 

where 7 -I is a positive constant. 

A traveling wave moving in x-direction with speed c is a 
solution of (1), (2) and (3) in which each of p, v and p can 
be written as a function of the waveform parameter y=x-ct. 
Taking g-0 and r= 0, the conservation law equivalents of ( 8 ), 


(9) and (10) are 


<L£l 

at 


.fl- 

ax 


(pv) 


fl_ 

3 t 


(PV) = -IT(PV 2 + P - «f“> 


£v. 


fl- 

at 


( hpv 2 +pu) 


ax 

a_ 

ax 


ax' 

[ h p V 2 + p + P U ) V ] + J 7 ( V 7J ) 


a r 3 . 3 v. . 


(ID 

( 12 ) 

( 13 ) 


Denoting the conditions before and after shock with 
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subscripts 1 and 2 respectively, (11), (12), (13) yield 


P 1 (V X -c) - p 2 (v 2 -c) (14) 

Pl(Vi-c) 2 + Pi - - P 2 (v 2 -c) 2 + P 2 - V 2 fa Z (15) 

l%Pl (Vj - 0)2+?!+^^! } (V X - c) - ( V J - c )»?i§7 1 -' c §7('?lf7 1 ) “ 

{Hp 2 (v 2 - C ) 2 +p 2 +p 2 u 2 ) (v 2 -c) - (v 2 - c )»?2fj 2 -' c f7(»?2fj 2 > ( 16 ) 


As mentioned before, a traveling wave moving in x- 
direction with speed c is a solution of (1), (2), (3) in 
which each of p , v and p can be written as a function of the 
waveform parameter y-x-ct . The waveforms of such a wave 


differ from those derived from conventional theory. Hence 
an alternative to the Rankine - Hugonio t equation [4, p.317] 


may be desirable. From (14), (15), (16) one can show that 

(u 2 +p 2 p 2 ) - (u 1 +p 1 p 1 )-*S(p 2 - Pi ) (p 2 1 +Pl 1 )+h(p ~ 2 1 -Pi 1 ) 

. ( „ iv 2 . K £v 2 3—, 3-¥-l \ -I / 1 7 \ 

(,2 ax ni dx ) + ^J^3x (,2 dx ^'dx^dx ^ (17) 

When ^ 2 "^!“® and /c-0 , (17) reduces to the conventional 


equation . 


4. Proposed Continuation of this Research 

The assumption of pressure diffusion leads to seemingly 
realistic waveforms for one dimensional traveling waves with 
shocks. However, an accompanying mechanism for calculating 
waveforms through a shock remains to be devised. 

If, as indicated above, a shock is very cold, metastable 
phase transition phenomena might be involved and a local 
alternative specification of internal energy u(p,p) must be 
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made. A first step is the use of the Van der Waals gas 
model. The internal energy differential for a Van der Waals 
gas is 


du 


C < - b/>) 


(?-•> 


1 - Z 


1 

(3 + a^ 1 - J 





(18) 


The proposed continuation of this research involves 
use of the pressure diffusion theory to devise a means of 
calculating p, v, and p values through a shock ( conserving 
mass, momentum and energy) starting with du as given by (18). 
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